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Abstract. We present the azimuthally averaged radial 
distribution of 748 regions of OB star formation in the 
whole galactic disk, based on our previous CS(2-1) sur- 
vey of UC H II regions. Embedded massive stars produce 
a total FIR luminosity of 1.39 10 s Lq within the range 
0.2 < R/R < 2 in galactocentric radius. We hnd 492 mas- 
sive star forming regions within the solar circle, producing 
81% of the total FIR luminosity. Separate analyses of the 
349 sources in the I and II quadrant (north), and of the 
399 sources in the III and IV quadrant (south), yield FIR 
luminosities (extrapolated to the complete galactic disk) 
of 1.1710 8 L Q and of 1.60 10 s Lq, respectively. Massive 
star formation is distributed in a layer with its centroid 
Z (R) following that of molecular gas for all galactocentric 
radii, both north and south. Its thickness for R < R a is 
~ 73 pc (FWHM), 62% the thickness of the molecular gas 
disk. The FIR luminosity produced by massive stars has 
a well defined maximum at R = 0.55 R , with a gaussian 
FWHM of 0.28 R - compared with 0.51 R Q for the H 2 sur- 
face density distribution. Toward the outer Galaxy, down 
from the maximum, the face-on FIR surface luminosity 
decays exponentially with a scale length of 0.21 R a , com- 
pared with 0.34 R for the H 2 surface density. Massive star 
formation per unit H 2 mass is maximum for R ~ 0.55 R 
in the southern Galaxy, with a FIR surface luminosity to 
H 2 surface density ratio of ~ 0.41 Lq/Mq, compared with 
~ 0.21 Lq/Mq at the same radius in the north, and with 
an average of ~ O.I8L0/M0 for the whole galactic disk 
within the solar circle. 
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medium: H II regions — interstellar medium: molecules, 
millimeter lines radio lines: molecular 



Send offprint requests to: L. Bronfman 

* Based partly on results collected at the European Southern 
Observatory, La Silla, Chile 



1. Introduction 

How many regions of massive star formation are there in 
the galactic disk and how are they distributed? Do all 
these regions have similar luminosities? Dust extinction 
prevents optical observations of OB associations farther 
than 3 kpc from the Sun. An alternative way to study the 
population of recently formed OB stars in the Galaxy is 
to use observations at longer wavelenghts in the far in- 
frared (FIR), where the ISM is more transparent. Ultra- 
violet light from OB stars, born in the dense dusty cores 
of giant molecular clouds, ionizes the closeby gas creating 
ultracompact (UC) H II regions and heating the surround- 
ing dust, which reradiates the energy mostly in the FIR. 
The star- forming region becomes thus detectable by IRAS, 
generally as a point-like source. But which of the many 
IRAS point-like sources in the galactic plane are heated 
by embedded massive stars? The FIR spectral characteris- 
tics of UC H II regions identified with the VLA have been 
studied by Wood & Churchwell (1989a), who calibrated 
their FIR colours as a tool for their identification. Apply- 
ing their criterion to the IRAS point-like sources in the 
galactic plane, they found 1646 sources which, according 
to their IRAS 60/i/12/x and 25/i/12/i flux ratios, are can- 
didates to be sites of OB star formation containing UC H 
II regions (Wood & Churchwell 1989b; WC89). 

Millimeter-wave emission lines from dense gas cores, 
where OB stars form, can be used to obtain their kine- 
matic distances and derive the FIR luminosity of their as- 
sociated dust, heated by the embedded stellar sources. The 
12 CO(l — 0) line has been surveyed toward IRAS point- 
like sources in the outer Galaxy (Wouterloot & Brand 
1996; and references therein), using less restrictive FIR 
colors criteria. Within the solar circle, however, where the 
majority of massive star formation takes place, there are 
several molecular clouds at different distances for most 
lines of sight, yielding complex CO velocity profiles. The 
12 CO(l — 0) line, affected by large optical depths, will tell 
us little about the dense interior of molecular clouds. The 
13 CO(1-0) and C 18 O(l-0) lines, of intermediate and low 
opacity in the galactic plane, are still sensitive to molec- 
ular gas column density, which can be enhanced not only 
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Fig. 1. Distribution of massive star forming regions in the whole Galaxy: (a) Longitude- velocity; and (b) Longitude- latitude. 
Sources with FIR fluxes higher than 3 10~ 10 Wm~ 2 , corresponding to an 08 star with a luminosity of 1.7 10 5 Lq at a distance 
of 4.25 kpc, are represented with open circles 



toward dense cores but also toward extended regions of 
intermediate density. 

One of the best specific tracers of dense gas is the CS 
molecule. Its J = (2 — 1) rotational line, at a frequency of 
99 GHz, has an excitation density threshold of 10 4 - 10 5 
cm~ 3 and is ubiquitous toward regions of massive star for- 
mation but seldom elsewhere in the galactic disk - with the 
exception of the galactic center region (Bally et al. 1987). 
A complete survey of the CS(2 — 1) line toward IRAS 
point-like sources with FIR colors of UC H II regions in 
the galactic plane (Fig. 1) was presented by Bronfman et 
al. (1996; BNM). Nearly in all cases the CS(2 - 1) profiles 
obtained in such survey are singular; the dense gas cores 
are associated one-to-one with the heated dust which pro- 
duces the FIR emission detected by IRAS as a point-like 
source. The CS(2 — 1) line, as well as other higher fre- 
quency rotational lines of CS, have been used successfully 
to map dense molecular gas cores and to derive their physi- 
cal parameters (Plume et al. 1997; and references therein). 

The main goal of the present paper is to derive the 
mean radial distribution of the facc-on FIR surface lumi- 



nosity produced by embedded massive stars in the whole 
Galaxy, and to compare it with the mean radial distribu- 
tion of H2 surface density as derived from CO surveys of 
the Milky Way. For a direct comparison, thus, we analyze 
our sample of OB star formation regions in a manner as 
close as possible as that used for the derivation of the H2 
surface density (Bronfman et al. 1988; Paper I). 

As an initial step in the analysis of the BNM CS(2 — 1) 
survey, we present the azimuthally averaged distribution 
of regions of massive star formation along the whole galac- 
tic disk, 0.2 < R/R < 2. We also perform separate anal- 
yses of the northern (I and II quadrants) and southern 
(III and IV quadrants) Galaxy in order to search for large 
scale deviations from axial symmetry. Using the kinematic 
distances derived from the CS(2 — 1) line profiles and the 
FIR fluxes from the IRAS Point Source Catalog (PSC), 
we evaluate the FIR face-on surface luminosity, averaged 
over galactocentric rings 0.1 R (0.85 kpc) wide, for the 
northern, southern, and complete galactic disk. The re- 
sults obtained are compared with the molecular hydrogen 
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Fig. 2. Distribution of detected and undetected sources integrated in (a) galactic latitude and (b) galactic longitude 



surface density distribution as derived from large scale 
CO(l — 0) surveys of the Galaxy. 

In Sect. 2 we summarize the observations and the re- 
sults of the BNM CS(2 - 1) survey. In Sect. 3 we de- 
scribe our axisymmctric analysis and derive the mean ra- 
dial distribution of OB star formation in the Galaxy, the 
azimuthally averaged face-on FIR surface luminosity these 
stars produce, and the mean FIR luminosity per OB star 
forming region. We evaluate the shortcomings in our anal- 
ysis which are due to the spiral structure of the Galaxy and 
to the incompleteness of our sample. In Sect. 4 we compare 
these results with the surface density of cold molecular gas 
as traced by CO, discuss the large scale asymmetries of OB 
star formation in the Galaxy, and estimate the contribu- 
tion of embedded OB stars to the FIR luminosity of Giant 
Molecular Clouds (GMCs). In Sect. 5 we summarize our 
conclusions. 

2. Observational results 

The BNM CS(2 - 1) survey was carried out between 1989 
and 1992 with the 15-m SEST (Swedish-ESO Submillimc- 
ter Telescope) at La Silla, Chile, and with the OSO (On- 
sala Space Observatory) 20-m Telescope, Sweden, in May 
1993. The angular resolution of the telescopes at 99 GHz, 
the CS(2-1) line frequency, is 50" for SEST and 39" for 
OSO (main beam FWHM), similar to the IRAS survey 
position accuracy. The instruments, observing techniques, 



and data calibration are described by BNM. The observed 
sources were selected from the IRAS PSC (Version 1), 
using exactly the same criterion as in WC89. They are 
all located near the galactic plane, they have a 25/x/12/x 
flux ratio larger than 3.7 and a 60/i/25/i flux ratio larger 
than 19.3, and are not identified as stars or galaxies. A 
total of 1427 sources were selected. The survey covered 
—2° < b < +2° for longitudes corresponding to the in- 
ner Galaxy, 300° < I < 360° and 0° < I < 60°. For 
the rest of the galactic plane the latitude coverage was 
-4° < b < +4°. 

Of the 1427 sources observed, the survey detected 843 
(59%) above a level of - 0.3 K peak T A (BNM), with 
a few detections below this level, particularly in the far 
outer Galaxy (Table 1). Of the total number of sources 
detected, 37% arc in the I quadrant; 10% in the II quad- 
rant; 11% in the III quadrant, and 42% in the IV quad- 
rant. The fraction of detected sources is 0.58 toward the 
outer Galaxy (II and III quadrants), and 0.60 toward the 
inner Galaxy (I and IV quadrants); toward certain longi- 
tudes in the direction of the molecular annulus, however, 
the number of detected sources more than doubles that of 
undetected ones. 

The distribution in longitude of the detected sources, 
integrated over latitude, shows two prominent maxima 
around / = 30° and I = 330° , toward lines of sight tangent 
to a galactocentric circle 0.5 R in radius (Fig. 2a). These 
maxima, symmetric with respect to I = 0°, arc strong 
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Table 1. Distribution in galactocentric radius of the 
sources included in the analysis 
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(25) 
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(15) 
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(18) 
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(24) 
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1.15 
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(24) 


23.30 


1.25 
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16 
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25.38 


1.35 
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(22) 
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27.46 


1.45 
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29.54 


1.55 
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31.62 


1.65 
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33.70 


1.75 
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1 
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1.85 


1 


(1) 


2 


(2) 


37.87 


1.95 


1 
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(0) 


39.95 


Total 


349 


(272) 


399 


(304) 


394.1 



a In parenthesis the number of sources when a lower lumi- 
nosity cutoff is applied to the sample (see Sect. 3.5) 



direct observational evidence for the existence of an an- 
nulus of massive star formation in the Milky Way that is 
coincident with the molecular annulus, as traced by CO. 
The center of the distribution in galactic latitude of the 
detected sources (Fig. 2b) is very close to b — 0°, with 
an average of -0?05 ± 0?03 and a FWHM of - 0?8. The 
IRAS point-like sources not detected in CS, in compari- 
son, are more evenly distributed in the galactic disk, and 
their latitude distribution is about 20% wider. 

3. Mean radial distribution of embedded OB stars 

Using the kinematic information provided by the CS(2 — 1) 
line for the detected IRAS sources, we have derived their 
galactocentric distances assuming purely circular motion 
about the galactic center. For simplicity we used a linear 
rotation curve and the standard IAU constants, R — 8.5 
kpc and Vq = 220 km s _1 . For lines of sight toward the 
galactic center and anticenter the VLSRs of the sources 
are small and may have large fractional errors due to sys- 
tematic or random peculiar motions. Thus, for the anal- 
ysis here we have excluded all sources within 10° of the 
galactic center and within 5° of the anticenter. Therefore, 
a total number of 748 detected sources have been used for 
the present analysis (Table 1), of which 349 (47%) are in 
the northern Galaxy (I and II quadrants) and 399 (53%) 
arc in the southern Galaxy (III and IV quadrant). 



After finding the galactocentric distance R for each 
source, they were binned in annuli 0.1 R wide. The total 
number of sources in each ring was divided by the ring 
sampled area to evaluate the azimuthally averaged num- 
ber surface density of massive star formation regions as 
a function of galactocentric radius north, south, and for 
the complete dataset (Table 2; Fig. 3). The derivation de- 
pends only on the assumption of pure circular motion, and 
uses the kinematic galactocentric distance of each source, 
which can be derived one-to-one from the CS(2-1) velocity 
profiles. 

3.1. Centroid and thickness of the OB star formation 
layer 

To derive the azimuthally averaged centroid Z (R) and 
thickness Zi (i?) (HWHM) of the massive star formation 
layer we need to calculate the first and second order mo- 
ments of the Z distribution for each of the galactocentric 
rings analyzed. Within the solar circle, however, for each 
source there are two possible kinematic heliocentric dis- 
tances allowed by the rotation curve and, therefore, two 
possible values of Z. Following Paper I, we assume that 
the distribution P(Z) of sources is of gaussian form in Z, 

,r*s , \n2\Z- Z (R)} 2 
P(^)=Aexp (!) 

where R is the galactocentric distance, Z D (R) and Zi(R) 
are fixed for each ring and A is a normalization constant. 
We then calculate the first and second order moments of 
the Z distribution in each ring, Zi(R) and Z 2 {R), by eval- 
uating a statistical average, using weights that depend on 
the near or far distance from the plane, Z n and Zf. 

N(R) 

Zi(R) = £ [P n (j) Z n (j) + P f (j) Z f {j)] (2) 

and 

N{R) 

Z2{Rf = Wr~)^ [PnU) ZnU)2 + Pf(j) Zf{j)2] (3) 
j=i 

where N(R) is the number of sources in a ring at galac- 
tocentric distance R, P n = P(Z n ), Pf = P(Zf), and the 
normalization condition P n + Pf = 1 is used. 

The centroid Z (R) is equal to the first order moment 
Z^R), and the thickness Zi(R) (HWHM) is given by 

Zi (R) 2 = 2 ln2 [Z 2 (R) 2 - Z (R) 2 } (4) 

The weights P n and Pf depend in turn of Z Q and Zi 
at each ring, so the problem is cyclic. It can be solved 
by an iterative procedure (Casassus 1995), starting with 
reasonable values, like Z„ = and Zi = 60 pc, close to 
the average HWHM of the H 2 layer within the solar circle. 
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Table 2. The embedded massive stars layer 
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154 ±6 
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0.63 ±0.16 
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-84 ±4 
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94 ±9 


127 ±9 


1.35 
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133 ±8 


157 ±5 
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71 ±5 


265 ±6 
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While it is clear that each source is either at the near 
or at the far distance within the solar circle, the weights 
used in the statistical averages account for the probability 
of each source to be at the near or far distance by using 
the well known latitude effect, i. e., the farther away a 
source, the more likely to be close to the galactic plane. 
The weights are used only to derive large scale averages in 
a statistical sense, and not individually to derive a mean 
distance for each source. The results of the iterative calcu- 
lation are displayed in Fig. 4 and listed in Table 2; empty 
spaces in the table correspond to those rings where the 
number of detected sources is less than 3, not enough to 
have a convergent iteration so as to derive values for Z Q 
and Z i . 

2 

3.2. Azimuthally averaged FIR surface luminosity 
produced by embedded OB stars 

The azimuthally averaged face-on FIR surface luminosity 
is evaluated by summing the luminosities of the sources in 
each ring and dividing the result by the sampled area. The 
average luminosity per source is the ratio of the summed 
luminosity in each ring to the number of sources consid- 
ered. We estimate the FIR flux of a source as the sum over 
the four IRAS bands, 

4 

Firas = ^F v {j), (5) 

where F v (j) are the IRAS band flux densities, as listed 
in the IRAS PSC (Pcrault 1987). The FIR luminosity for 
each source, outside the solar circle, is computed as L = 



4tt D 2 Fjras , where D is its kinematic distance, obtained 
from its CS(2-1) line profile. We overcome the two-fold 
distance ambiguity, within the solar circle, by evaluating 
the total luminosity in each ring, L(R), as a weighted sum. 
Like for the Z moments, in that sum we use the latitude 
effect to take into account the probability for each source 
to be at the near or far distance in the galactic disk, i. e., 

N(R) 

L(R) = l p nti) L n (j) + Pf(j) L f (j)} (6) 

where L n = 4tt Firas and Lf = An D"j Fjras are the 
luminosities that the source would have at the near and far 
distances, respectively. The weights P n = P(Z n ) and Pf = 
P(Zf) are obtained using the values of Z D and Zi listed in 
Table 2. For the innermost ring, at R/R = 0.25, where no 
model is available, and at Rj R — 0.35 in the north, we use 
the subcentral distances to estimate the luminosities. The 
FIR surface luminosity produced by regions of massive 
star formation and the average FIR luminosity per region, 
as a function of R, are shown in Table 3 and Fig. 3. The 
procedure used to remove the distance ambiguity, within 
the solar circle, is fully equivalent to the direct partitioning 
technique used to derive the H 2 surface density in Paper 
I (see Eq. (4) and Fig. 15). 

3.3. Errors derivation 

The errors in the number surface density consider only 
the counting statistics assuming a Poisson distribution. 
The uncertainties in the best values of the parameters Z a 
and Zi were estimated on the basis that the errors are 

2 
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Fig. 3. The left column shows the number of detected sources divided by the sampled area, as a function of R, 
for the northern (0° < I < 180°) Galaxy (top row), the southern (180° < I < 360°) Galaxy (medium row), and 
the whole Galaxy (bottom row). The center column shows, in the same order, the azimuthally averaged face-on FIR 
galactic surface luminosity produced by embedded massive stars; the right column shows the azimuthally averaged 
FIR luminosity per embedded massive star forming region. The solid line corresponds to the complete dataset. The 
dashed line represents the dataset corrected for incompleteness of the sample, as described in Sect. 3.5 



dominated by the peculiar motions of the sources, which 
we assume to be of ~ 5 km s _1 (VLSR). Using this value 
we have then calculated the uncertainties in the kinematic 
distances, and in the quantities which are derived from the 



kinematic distances. At 1.1 < R/R < 1.2, in the northern 
Galaxy, the average luminosity per source is dominated by 
one source, G70.293 (BNM), a very bright UC H II region 
in the Perseus spiral arm. 
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Fig. 4. The azimuthally averaged centroid Z (top panel) and the thickness (HWHM) Zi of the massive star formation 
layer (bottom panel) arc shown for the northern (0° < I < 180°) Galaxy (left); the southern (180° < I < 360°) Galaxy 
(center); and the complete Galaxy (right) 



Table 3. FIR surface luminosity 
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The uncertainties in the Lpjn values have also a con- 
tribution from the instrumental errors in the IRAS detec- 



tors. We used the errors in the 100 /jm band as repre- 
sentative of the fractional flux uncertainty for each source 
(Casassus et al. 2000). When the flux value obtained by 
IRAS is only an upper limit, the error has been assumed 
to be of —100%. There is also some level of correlation 
in the position of the sources; as an example, for mas- 
sive star forming GMCs in the IV galactic quadrant there 
are ~ 4 embedded UC H II regions per GMC (Bronfman 
1992). We have taken into account such correlation by 
computing for each galactocentric ring N corr (R), the av- 
erage number of sources that lie within a box of sides 2° 
by 2° by 20 km s -1 - a representative maximum value for 
a GMC - and multiplying the errors in Lpjji at each ring 
by ^jN corr (R). 

3.4- Spiral structure and the axisymmetric analysis 

What are the shortcomings in our analysis caused by devi- 
ations from axial symmetry in the galactic distribution of 
embedded OB stars? Images of external galaxies clearly 
show a high concentration of massive star formation in 
spiral arms, and the same should be expected for our own 
Galaxy. While the large scale distribution of OB stars 
in the Milky Way cannot be traced optically because of 
dust obscuration in the disk, H II regions produced by 
these stars can be observed in radio-recombination lines, 
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Velocity LSR [Km/s] 

Fig. 5. The spiral arm pattern used for the analysis. The 
arms are labelled 1 (Carina), 2 (Crux), 3 (Norma), and 
4 (3-kpc). The CO emission (Paper I; Grabelsky et al. 
1987) integrated in latitude from —2° to 2° is represented 
in gray-scale. Lowest contour level is at 3 a. The scale in 
the inset figure is in kpc, and the galactic center is labelled 
G.C. 



like H 109a, and have been used to study the spiral arm 
structure of our Galaxy (Georgelin & Georgelin 1976; 
Downes et al. 1980). The UC H II regions, which point at 
very young star forming regions, still embedded in their 
parental molecular clouds, could be even more confined 
to spiral arms than the extended H II regions observed in 
centimetric wavelengths. 

Spiral arms within the solar circle appear to be more 
separated and less difficult to trace in the southern Galaxy 
than in its northern counterpart. These arms, as outlined 
by CO emission from molecular clouds, are tangent to the 
line of sight at longitudes 283° (Carina); 308° (Crux); 328° 
(Norma); and 337° (3-kpc expanding arm). The tangent 
longitudes are traced by discontinuities in (a) the galactic 
CO emission integrated in velocity and latitude 1(1) (Pa- 
per I; Grabelsky et al. 1987) and in (b) the rotation curve 
derived from CO observations (Alvarez et al. 1990). The 



distribution of H II regions in the southern Galaxy, de- 
rived from an extensive survey of hydrogen recombination 
lines toward 5 GHz continuum sources (Caswell & Haynes 
1987), is consistent also with the presence of these four 
spiral arm segments in the IV galactic quadrant. 

To test our axisymmetric analysis we have derived the 
mean radial distribution of UC H II regions in the south- 
ern Galaxy in the framework of this generally accepted 
spiral arm model. Most GMCs identified through CO ob- 
servations in the southern Galaxy are coincident, in space 
and velocity, with one or more radio H II regions from the 
Caswell & Haynes (1987) sample, as well as with several 
UC H II regions from the dataset we analyze here. We 
have used those radio H II regions from Caswell & Haynes 
that have well denned distances, derived through H 2 CO 
absorption lines, to remove the distance ambiguity for the 
associated GMCs and for their embedded UC H II regions 
(Bronfman 1992). These GMCs, together with the tan- 
gent points listed above, trace a spiral arm pattern in the 
longitude- velocity diagram (Fig. 5). The distance ambigu- 
ity for the remaining GMCs is removed by assigning each 
one to the closest spiral feature in the longitude-velocity 
diagram. For some GMCs in the far side of the Crux arm, 
close to the solar circle, we use their size-to-linewidth ratio 
to distinguish them from local clouds with small VLSRs. 

The mean radial distribution of FIR surface luminosity 
produced by embedded UC H II regions, obtained using 
the spiral model, is fairly similar to that obtained using the 
axisymmetric analysis (Fig. 6a). The average luminosity 
per source (Fig. 6b) shows a secondary peak, close to R = 
R , caused by the far side of the Crux spiral arm cutting 
through the solar circle at 316° < I < 320°. We do not 
quote errors in this comparison because the spiral model, 
although well updated, involves much uncertainty and we 
are mainly interested in the general trend of the mean 
radial variations. It is apparent that the peaks are at the 
same positions and the general shapes are the same for 
both analyses. The axisymmetric analysis, however, seems 
to overestimate the mean FIR luminosities for radii other 
than the main peak locus. 

3.5. Completeness of the sample and contamination by 
low luminosity sources in the solar neighborhood 

The question of completeness has been adressed in general 
by WC89, who show that UC H II regions are tightly 
confined in the FIR color-color plot, and can be easily 
distinguished from other entries in the IRAS PSC. Thus, 
the distinctive FIR colors of embedded massive stars can 
be used to count the number of such objects in the Galaxy. 
They conclude that all embedded O stars and some B stars 
will be detected by the color-color criterion applied to the 
IRAS PSC in the Galaxy. Contamination of the sample 
by cloud cores with lower mass stars has been analyzed 
by Ramesh & Sridharan (1997), who argue that the total 
number of potential UC H II regions in the WC89 sample 
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Fig. 6. (a) FIR surface luminosity and (b) average lu- 
minosity per source, as a function of R, for the southern 
Galaxy. The solid line shows the results obtained by re- 
moving the two-fold distance ambiguity, within the solar 
circle, using a spiral model. The dotted line shows the 
results obtained using the axisymmetric analysis 

is only about 25% of the sample. Our own criterion, based 
on the detection of the CS(2-1) line toward each of the UC 
H II region candidates (BNM), brings down the original 
WC89 sample by a factor of <~ 2. 

We analyze here the question of completeness of the 
sample as a function of distance; this is particularly nec- 
essary to understand the behaviour of the number density 
and mean luminosity of sources near the solar circle. The 
FIR luminosity versus distance for all sources in the sam- 
ple analyzed is shown in Fig. 7. The fairly homogeneous 



2 4 6 
D/kpc 




Fig. 7. FIR luminosity versus distance for all the sources 
in the sample. Within the solar circle (left) the dots rep- 
resent distances and luminosities obtained through the 
weighting method (See Sect. 3), while the circles repre- 
sent subcentral sources, with better determined distances. 
Outside the solar circle (right) all distances are well de- 
fined 

distribution in distance of sources with FIR luminosities 
larger than log(L/L Q ) ~ 4 suggests that our sample is 
complete above such luminosity level. Closer than ~ 3 kpc 
from the Sun the sample appears to be contaminated by 
lower luminosity sources. To test our results we have rean- 
alyzed our data applying a lower FIR luminosity cutoff to 
the sample, both for the spiral and for the axisymmetric 
model. Wouterlout et al. (1995, their Fig. 20b) show that 
the WC89 FIR color criterion for UC H II regions tends 
to select point sources with log(L/L@) > 4. Thus, even if 
the IRAS PSC were more sensitive, there would be little 
point in reducing the luminosity limit much under 10 L@. 

The CS(2-1) detection requirement appears not to in- 
troduce additional biases so we consider unnecessary to 
address separately the completeness of the CS survey. The 
observed antenna temperatures Ta within the solar circle, 
in particular for the subcentral sources, are not distance 
dependant. For sources outside the solar circle, the average 
Ta is 0.9 K for distances between 2 and 6 kpc, and 0.6 K 
beyond. The SEST bcamwidth of 50" (FWHM) projects 
a linear size of 2.4 pc at 10 kpc, about the typical size 
of CS cores harbouring massive star formation (Bronfman 
1992), so beam-dilution is not a stringent problem for the 
CS detections. 

The mean radial distribution of the FIR surface lu- 
minosity does not change noticeably when a lower cut- 
off luminosity is applied to the sample. We have com- 
puted the azimuthally averaged FIR surface luminosity 
and the mean luminosity per source as a function of galac- 
tocentric radius, using the same spiral analysis as in Sect. 
3.4 within the solar circle, and a lower cutoff luminosity 
log(L/L Q ) = 3.9. The results are presented in Fig. 8, and 
the number of sources above the cutoff is listed, for each 
ring, in Table 1 (in parenthesis). At the cutoff value used, 
the three peaks in the average luminosity per source have 
the same magnitude, i.e Lfir ~ 0.35 10 6 i Q . These peaks 
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Fig. 8. (a) FIR surface luminosity and (b) average lumi- 
nosity per source, as a function of (R), in the southern 
Galaxy. The two-fold distance ambiguity within the solar 
circle, has been removed using a spiral model. The results 
obtained using the complete dataset are shown with the 
solid line; results obtained when leaving out all sources 
with FIR luminosities lower than 8000 Lq are shown by 
the dotted line 

arc dominated by sources located in the far sides of the 
Norma arm (R = 0.55 R ), the Crux arm (R = 0.95 R ), 
and the Carina arm (R = 1.25 R ). 

While we cannot, for the time being, confirm that the 
average luminosity per source diminishes with galactocen- 
tric radius, it is tempting to suggest that UC H II regions 
appear to be very good tracers of spiral arm structure in 
the Galaxy, and that a full spiral analysis of their distribu- 



tion in the Milky Way is mandatory. While in the present 
paper we perform a zero'th order axisymmetric compari- 
son of the distribution of OB star formation regions and 
of molecular clouds, there is an evident necessity of resolv- 
ing the two-fold distance ambiguity for all these regions 
in order to refine the present results, within the solar cir- 
cle, and extend them toward a more exact analysis of the 
mean luminosity per source, which has a direct impact 
on the determination of the IMF for massive star forming 
regions. 

The FIR surface luminosity derived from the axisym- 
metric analysis, our main result, remains also fairly un- 
changed (Fig. 3) when a lower luminosity cutoff is applied 
to the sample. The FIR surface luminosity at each radius 
is the product of the number of sources considered, per 
unit area, times their average luminosity. While the num- 
ber surface density of sources decreases when a luminosity 
cutoff is applied, particularly for the solar neighborhood, 
the average luminosity per source increases, since we are 
not taking into account the less luminous sources for the 
average. Therefore the results we present next, based on 
our determination of the FIR surface luminosity produced 
by embedded massive stars and its comparison with the 
H 2 surface density within the solar circle, appear to be 
not too affected by problems of completeness. As for the 
effects of spiral structure over the derived FIR surface den- 
sity, they are small and would have presummably the same 
effects over the distribution of molecular hydrogen derived 
using an axisymmetric model. 

4. The galactic disk traced by embedded OB 
stars and by molecular clouds 

4-1. Centroid and thickness of the disk 

The formation of OB stars in the galactic disk is dis- 
tributed in a layer coincident but thinner than that of 
molecular clouds as traced by CO. Fig. 9 shows the radial 
variations of the centroid Z a and the thickness (HWHM) 
Zi of the embedded massive star layer for the southern 
and northern Galaxy, compared with their H2 counter- 
parts. The molecular gas parameters were derived from 
the Columbia U. - U. Chile deep CO surveys of the I and 
IV galactic quadrants ( Paper I; Cohen et al. 1986); from 
surveys of the Carina region by Grabclsky et al. (1987) and 
of the III galactic quadrant by May et al. (1988, 1993), 
analyzed by Casassus (1995); and from a composite CO 
dataset of the II quadrant analyzed by Digcl (1991). 

The average thickness of the embedded massive star 
distribution, within the solar circle, is 73 pc (FWHM), 
compared with 118 pc (FWHM) for the layer of molecular 
gas traced by CO (Tables 6 and 7). It has been shown by 
Dame & Thaddeus (1994), however, that there are pos- 
sibly two components for the molecular gas layer: a thin 
one, with a FWHM of 80 pc, similar to the thickness of 
the young massive stars layer derived here, and a thick 
and diffuse one, with a FWHM in excess of 500 pc. 
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Fig. 9. The azimuthally averaged thickness (HWHM) Z± (top) and the centroid Z Q (bottom) of the massive star 
formation layer (solid line), as a function of R, are compared to those of H 2 (dashed line), as derived from the Columbia 
U. - U. Chile deep CO surveys of the Galaxy 



The warping of the disk outside the solar circle, as 
traced by the centroid Z Q of both the massive star forma- 
tion and the H2 layer, appears to be more pronounced in 
the northern Galaxy than in the south. A tentative expla- 
nation is that the spiral arms outside the solar circle can 
be traced within a larger longitude range in the south than 
in the north; we may be therefore averaging out the warp 
in azimuth in the south because of the large longitude 
span (Casassus 1995; May et al. 1997). For example, it is 
apparent from Fig. 2 than there are many more sources 
detected outside the solar circle in the IV quadrant than 
in the I quadrant of the Galaxy. 

The flaring of the molecular gas layer in the outer 
Galaxy is evidenced by both the molecular and by the em- 
bedded massive star distribution in the south. However, 
while the thickness of the massive star layer in the north- 
ern outer Galaxy grows like that in the south, the mean 
thickness of the molecular hydrogen layer in the north- 
ern outer Galaxy is only about half that of its southern 
counterpart. A possible explanation is that while the lon- 
gitude coverage of our survey of embedded massive stars 
is equally complete north and south, the composite CO 
data set of the northern outer Galaxy analyzed by Di- 
gel (1991) only covers the region 65° < I < 116°, com- 
pared to the complete longitude coverage of the III and 
IV quadrant CO surveys used by Casassus (1995) in his 



analysis. Therefore the azimuthal range for the average is 
much larger in the south than in the north; averaging the 
galactic warp over large scales would increase the mean 
thickness derived for the molecular gas layer. 

4-2. Embedded OB stars, FIR surface luminosity, and H 2 
surface density 

The mean radial distribution for the 748 massive star 
forming regions analyzed here (Fig 3) is maximum at 
R ~ 0.55 i? , falls sharply toward the galactic center, and 
decays exponentially toward the outer Galaxy (Fig. 10). 
Of these regions, 492 (67%) are found at R < R , and 
332 (44%) at 0.4 < R/R Q < 0.8, the region of the molec- 
ular annulus (Paper I). The fraction of 44% grows to 55% 
when the contamination of the sample by low luminos- 
ity regions in the solar neighborhood is considered (Sect. 
3.5). A similar annular radial dependence for the number 
surface density of OB associations in the inner Galaxy 
has been found by Mc Kee & Williams (1997) using an 
older compilation of radio H II regions (Smith et al. 1978) 
binned in rings 2 kpc wide. Such annular shape was found 
also by Kent et al. (1994) for the NIR emission attributed 
to M giants, based on large scale satellite observations. 

The FIR surface luminosity produced by embedded 
massive stars for the Galaxy as a whole is even more re- 
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Fig. 10. From top to bottom: azimuthally averaged face- 
on FIR surface luminosity, as a function of R, originated 
by embedded massive stars (solid line) for the northern 
Galaxy (0° < I < 180°); the southern Galaxy (180° < I < 
360°); and the whole Galaxy, compared with molecular 
hydrogen face-on surface density (short dashed line), as 
derived from the Columbia U. - U. Chile deep CO surveys 
of the Galaxy. Data are uncorrected for the expansion ve- 
locity field at the 3-kpc arm locus (Paper I) 



stricted to the molecular annulus than their number den- 
sity. Integration over the whole galactic disk yields a total 
FIR luminosity of 1.39 10 s L . Within the solar circle the 
FIR luminosity is 1.131O 8 L , 81% of the total, and for 
0.4 < R/R < 0.8, the FIR luminosity is 0.97 10 s L Q , rep- 
resenting 70% of the total. The radial dependance of the 
FIR surface luminosity produced by embedded massive 
stars has also a narrower maximum than the H2 surface 
density. A gaussian fit to the azimuthally averaged FIR 
surface luminosity, within the solar circle, yields a FWHM 
of 0.28 i? (Table 4), - 55% the FWHM obtained for the 
H 2 radial distribution, 0.51 R (Table 5). Down from the 
maximum, toward the outer Galaxy, the FIR surface lu- 
minosity decreases, as a function of R, much faster than 
the H 2 surface density, with an exponential scale length 
of 0.21 R , compared to 0.34 R for the H 2 . 

The ratio of FIR luminosity produced by embedded 
OB stars to H 2 mass, which may be used as a large 
scale tracer of massive star formation efficiency, is not 
constant throughout the Galaxy and has also a peak at 
R ~ 0.55 R . Integrating the H 2 surface density, as traced 
by the Columbia U. - U. de Chile CO surveys of the 
Galaxy, we compute a total H 2 mass of 6.11O 8 M for 
0.2 < R/R < 1, and of - 2 10 8 M Q for 1 < R/R Q < 
1.7. For all the CO data we have used main beam antenna 
temperatures T m b, with a main beam efficiency of 0.82 
(Paper I); CO to H 2 conversion factor A r (H 2 )/14 7 co = 
1.56 1 20 molecules cm" 2 (K km s" 1 )" 1 derived using the 
same CO data and EGRET gamma-ray observations by 
Hunter ct al. (1997); a distance of the Sun to the galac- 
tic center of R Q = 8.5 kpc; and no correction for helium 
abundance. The average FIR luminosity to H 2 mass ra- 
tio, within the solar circle, is ~ 0.18 L Q /M©; for 0.4 < 
R/Ro < 0.8 the ratio grows to ~ O.23L /M ; and for 
0.5 < R/Ro < 0.6 there is a maximum of ~ 0.29 L Q /M Q . 
Outside the solar circle, for R/R < 1.7, the derived mean 
ratio is ~ 0.13 L /M ; this figure can be considered as an 
upper limit, since the conversion factor iV(H 2 )/lFco has 
been argued to be from 2 to 4 times larger in the outer 
Galaxy (May et al. 1997). 

The radial distribution for the complete set of OB star 
formation candidate sites from WC89 has been derived by 
Comeron & Torra (1996), assuming axial symmetry but 
without the use of any kinematic information. Our present 
results do not confirm the peak in massive star formation 
that Comeron & Torra (1996) find at R = 2 kpc, although 
it is fair to say that our analysis is oriented to the galactic 
disk, and that we have therefore cut out the region within 
10° of the galactic center. However, we detect only 1 source 
in the range 1.70 kpc < R < 2.55 kpc, where we sampled 
50% of the face-on area, while for their result to hold there 
should be more than 50 OB star forming regions in the 
area we sampled. 
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Table 4. Azimuthally averaged face-on FIR surface luminosity produced by regions of massive star formation 



Sampled region 


Gaussian maximum 8 

(£q P c- 2 ) 


Gaussian center a 

(Ro) 


Gaussian FWHM a 

{Ro) 


Exponential scale-length b 

{Ro) 


North 


1.02 ±0.16 


0.60 ± 0.02 


0.31 ±0.06 


0.21 ± 0.03 


South 


1.59 ±0.16 


0.60 ± 0.01 


0.26 ±0.03 


0.22 ± 0.02 


Complete 


1.29 ±0.15 


0.60 ± 0.02 


0.28 ± 0.04 


0.21 ±0.02 



a Fit between 0.2 R a < R < R a 
b Fit between 0.5 Ro < R < 2 R a 



Table 5. Azimuthally averaged face-on surface density of molecular hydrogen 



Sampled region 


Gaussian maximum 8 

(M Q pc- 2 ) 


Gaussian center 8 

{Ro) 


Gaussian FWHM a 
{Ro) 


Exponential scale-length b 
{Ro) 


North 


5.5 ±0.4 


0.57 ±0.01 


0.37 ± 0.03 


0.29 ± 0.02 


South 


3.9 ±0.2 


0.55 ± 0.02 


0.71 ±0.08 


0.40 ± 0.04 


Complete 


4.5 ±0.2 


0.57 ±0.01 


0.52 ±0.03 


0.34 ± 0.02 



a Fit between 0.2 R < R < R a 
b Fit between 0.5 R a < R < 2 R a 



4-3. Large scale asymmetries within the solar circle 

Among the most striking results of the early CO surveys 
of the Galaxy were the large scale north-south deviations 
from axial symmetry obtained from separate fits of ax- 
isymmetric models to the CO data (Robinson ct al. 1983; 
Cohen et al. 1985). The maximum molecular gas face-on 
surface density in the molecular annulus is a factor of 1.47 
higher in the I galactic quadrant than in the IV quad- 
rant (Paper I), where the maximum is really a plateau 
extending roughly for 0.4 < R/R < 0.8. Because massive 
star formation occurs mostly in giant molecular clouds, 
it would be expected to find a similar trend in the FIR 
surface luminosity produced by embedded massive stars. 
What we have found, however, is an enhancement of mas- 
sive star formation in the southern Galaxy within the solar 
circle. 

A fraction of 58% of the FIR luminosity produced by 
embedded massive stars, within the solar circle, comes 
from the southern Galaxy - while the total H2 mass, as 
traced by CO, is about the same north and south. A possi- 
ble explanation of such result could be that massive stars 
in the I galactic quadrant arc preferentially in the far side 
of the Galaxy, and/or that in the IV quadrant they are 
preferentially in the near side. Thus, by assuming axial 
symmetry in our models, we could be artificially inflating 
the luminosity in the south or lowering it in the north. 
However, because the FIR luminous dense molecular gas 
cores we have sampled in the CS(2-1) are well correlated in 
space and velocity with the population of giant molecular 
clouds in the Galaxy as traced by CO (Bronfman 1992), 
the same asymmetry should be found in the distribution 



of molecular gas. But what we find here for the FIR sur- 
face luminosity distribution is an asymmetry from north 
to south which is opposite in sign to that of the molecular 
gas surface density. 

Massive star formation per unit H2 mass, as a con- 
sequence, is higher in the IV than in the I galactic 
quadrant. The average FIR luminosity to H2 mass ra- 
tio in the southern Galaxy, for 0.2 < R/R < 1, is 
~ 0.21 L Cl) /Mq- for 0.4 < R/R < 0.8 grows to 0.28; 
and for 0.5 < R/R < 0.6 we find a galactic maxi- 
mum of ~ 0.41 Lq/M q . In comparison, the correspond- 
ing numbers in the northern Galaxy are ~ 0.16L Q /M Q , 
~ 0.18L©/M Q , and ~ O.21L /M . The enhancement 
of massive star formation at R ~ 0.55 R in the south- 
ern Galaxy appears to be correlated with the position of 
the Norma spiral arm tangent, where a large fraction of 
the most luminous sources in the sample are found. Mas- 
sive star formation enhancements of this kind have been 
observed in the spiral arms of external galaxies, like for 
instance in M51 (Nakai et al. 1994), although there the 
starburst seems to be driven by a companion galaxy. In 
our own galaxy such burst could be driven by a bar - that 
would have to be large enough to produce an effect at 
almost 5 kpc from the galactic center. 

4-4- The contribution of embedded OB stars to the FIR 
luminosity of GMCs in the Galaxy 

Preliminary results of our ongoing analysis of GMCs har- 
bouring embedded UC H II regions in the southern Galaxy 
show a fairly linear relation between the FIR luminosity 
from embedded OB stars and the virial mass of the clouds, 
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with a mean FIR luminosity to mass ratio of ~ 0.4 Lq /Mq 
(Bronfman et al. 2000). This ratio grows to <~ ALq/Mq, 
for GMCs associated with bright H II regions, when the 
FIR luminosity is computed from IRAS sky-flux images 
(Mooney & Solomon 1988). Therefore we suggest a lower 
limit of ~ 10% for the contribution of embedded OB stars 
to the total FIR emission from GMCs undergoing massive 
star formation. Our estimate can be considered as a lower 
limit because: (a) some embedded OB stars are likely to 
illuminate dust regions so large in angular size that they 
would not be identified as IRAS point sources and there- 
fore not included in our sample and, (b) if identified, some 
of the source FIR flux may have been lost because of the 
background subtraction. 

The ratio of FIR luminosity produced by embedded 
massive stars to virial mass for GMCs in the south- 
ern Galaxy, ~ 0.4 L Q /M Q , is similar to the FIR sur- 
face luminosity to H 2 surface density ratio found at the 
peak of the southern distribution of embedded OB stars, 
~ O.41L©/M . Therefore it is possible to suggest that 
most GMCs in that region of the Galaxy are presently 
forming massive stars, and it appears reasonable to use 
such value as a standard scale to evaluate the massive 
star formation efficiency per unit H 2 mass in the rest of 
the Galaxy relative to its maximum value. It is worth keep- 
ing in mind, though, that here we are averaging the H2 
surface density over large areas of the Galaxy, several kpc 2 
in size. Therefore our maximum value for the FIR lumi- 
nosity to H2 mass ratio will be much smaller than that 
obtained when computed only for the neighborhood of a 
massive star forming region, with a typical FIR luminos- 
ity of <~ 10 5 £© and a virial mass of ~ 410 3 M© (Plume 
et al. 1997) 

The total FIR luminosity of 1.13 10 8 L© derived here 
for the embedded massive star layer, for 0.2 < R/R Q < 
1.0, represents 6.2% of the FIR luminosity assigned to the 
H2 layer by Bloemen et al. (1990) from IRAS data, ~ 
1.83 10 9 Lq, and 5.5% of the FIR luminosity derived from 
COBE DIRBE data (~ 2.05 10 9 L©) by Sodrosky et al. 
(1997). Because in their respective analyses they consider 
all the FIR emission from the galactic plane, while here 
we compute only the FIR emission produced by dust fairly 
close to the heating stars, a strict lower limit of ~ 6% can 
be set for the contribution of embedded massive stars to 
the total FIR output from molecular gas within the solar 
circle. 



5. Conclusions 

Optical studies of massive stars are restricted to the so- 
lar neighborhood due to extinction by intervening dust. 
The FIR spectral signature of embedded OB stars pro- 
vides a unique opportunity to extend the study of their 
distribution to the whole galactic disk. We have analyzed 
here the first complete galactic survey of dense molecular 



Table 6. Separate axisymmetric models for the FIR lu- 
minosity from UC H II regions: R/R Q = 0.2 - 1.0 



Parameter 


North 


South 


Total FIR luminosity (£) 


9.43 10 7 Lq 


1.31 10 s Lq 


Mean surf. lum. < L > 


0.431/0 pc" 2 


O.6OL0 pc" 2 


Mean radius < R > 


0.64 R 


0.65 Ro 


Mean half- width < 2 1 > 


38 pc 


42 pc 


2 

Mean displacement < z > 


—6 pc 


—3 pc 



Table 7. Separate axisymmetric models for the H2 distri- 
bution as traced by CO emission: R/R = 0.2 - 1.0 



Parameter 


North 


South 


Total H 2 mass a (M) 


6.0 10 8 Mq 


6.2 10 s Mq 


Mean surface density < S > 


2.8M Q pc- 2 


2.9M pc~ 2 


Mean radius < R > 


0.65 Ro 


0.66 R 


Mean half- width < 21 > 


58 pc 


60 pc 


Mean displacement < z a > 


-12 pc 


3 pc 



a Using iV(H 2 )/Wco = 1.56 10 20 molecules cm" 2 (Kkms" 1 )" 1 
(Hunter et al. 1997) and R a = 8.5 kpc 



cloud cores associated with OB star formation. Our main 
conclusions are: 

1. We have derived the mean radial distribution for 748 
regions of massive star formation in the whole galac- 
tic disk. These regions produce a total FIR luminosity 
of 1.391O 8 L within the range 0.2 < R/R Q < 2.0. 
We find 492 regions, representing 67% of the sample, 
within the solar circle; they produce a FIR luminosity 
of 1.13 10 8 L Q , 81% of the galactic total. 

2. Separate analyses of the 349 sources in the northern 
Galaxy, and of the 399 sources in the southern Galaxy, 
yield total FIR luminosities (extrapolated to the com- 
plete galactic disk) of 1.17 10 8 Lq and of 1.60 10 s Lq, 
respectively. The fraction of the total FIR luminosity 
produced within the solar circle (~ 81%) is similar 
north and south. 

3. Regions of massive star formation in the galactic disk 
are distributed in a layer with its centroid Z (R) fol- 
lowing that of molecular gas for all galactocentric radii, 
both north and south. Within the solar circle the mean 
thickness of the OB star formation layer is ~ 73 pc 
(FWHM), a factor of 0.62 thinner than the H 2 layer 
(~ 118 pc FWHM). 

4. The FIR luminosity produced by massive stars, az- 
imuthally averaged over the whole galactic disk within 
the solar circle, has a well defined maximum at R — 
0.55 R , like the H 2 surface density, but with a gaus- 
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sian FWHM of 0.28 R a , a factor of 0.54 narrower than 
for the H 2 (0.52 R Q FWHM). 

5. Toward the outer Galaxy, down from the maximum, 
the face-on FIR surface luminosity, azimuthally aver- 
aged over the whole galactic disk, decays with an ex- 
ponential scale length of 0.21 R , steeper than the H 2 
surface density exponential decay which has a scale 
length of 0.34^. 

6. Massive star formation per unit H 2 mass is higher than 
average in the region of the molecular annulus. While 
the mean FIR luminosity to H2 mass ratio within the 
solar circle is - 0.18 L Q /M & , for 0.4 < R/R < 0.8 the 
ratio grows to ~ O.23L Q /M ; and for 0.5 < R/R < 
0.6 it reaches - 0.29 L Q /M Q 

7. Massive star formation per unit H 2 mass is maxi- 
mum in the southern Galaxy at 0.5 < R/R < 0.6. 
The average FIR luminosity to H 2 mass ratio there is 
~ 0.41 Z/ /M Q , compared with ~ 0.21 L Q /M Q for its 
northern counterpart. 

There is a fair amount of further work that we hope 
will follow from the present results. The availability of 
a fairly complete and homogeneous sample of OB star 
forming regions in the Galaxy allows statistical studies 
of the general physical process of massive star formation, 
and in particular of the high-mass end of the IMF. Bearing 
in mind that embedded massive stars should be the best 
tracers of spiral arm structure in our Galaxy, we foresee 
a concerted effort to resolve in a case-by-case basis their 
two-fold distance ambiguity within the solar circle. 
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